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A simple theoretical solution to the design of a control field that generates complete population 
transfer from an initial state, via N nondegenerate intermediate states, to one arbitrary member 
of M (M < N) degenerate states is constructed. The full control field exploits an (M + N — 1)- 
node null adiabatic state, created by designing the relative phases and amplitudes of the component 
fields that together make up the full field. The solution found is universal in the sense that it 
does not depend on the exact number of the unwanted degenerate states or their properties. The 
results obtained suggest that a class of multi-level quantum systems with degenerate states can be 
completely controllable, even under extremely strong constraints, e.g., never populating a Hilbert 
subspace that is only a few dimensions smaller than the whole Hilbert space. 

PACS numbers: 32.80.Qk 



I. INTRODUCTION 



Coherent control, i.e., controlling atomic and molecular processes via quantum interference effects induced by ex- 
ternal fields, has attracted great theoretical and experimental interest ^ n particular, our ability to manipulate 
atomic and molecular excitation and therefore drive the microscopic motion to generate an arbitrarily selected target 
state is one key prerequisite for quantum information processing jj] in atoms and molecules. 

One formal mathematical problem associated with coherent control is determining whether or not complete con- 
trollability is achievable. If, for a given system, it can be proved that there exists a control field (possibly composed 
of a number of component fields) that can yield complete population transfer to an arbitrary target state within a 
finite time, then the system is called completely controllable. For such systems, the primary task of coherent control 
studies is to design methods that are conceptually as simple as possible and also experimentally feasible. If the formal 
complete controllability of a quantum system is not established in general, then it is of great theoretical interest to 
construct specific control fields that can achieve, say, complete population transfer between two particular quantum 
states. 

Some existence theorems that establish conditions for complete controllability of a quantum system have been 
proved [H El 0- The Huang- Tarn-Clark theorem is believed to be the strongest result, but it applies only to 
systems with discrete and nondegenerate states. Ramakrishna et al. showed that for a quantum system with a 
Hilbert space of dimension L, the necessary and sufficient condition for complete controllability is that the field- free 
Hamiltonian and the interaction Hamiltonian induced by a control field generate a Lie algebra of dimension I? |(|. 
This criterion is applicable to multi- level systems that involve degenerate states |s, 9], but the required computations 
to generate the commutators of the Lie algebra structure can be demanding for large L, vary drastically from system 
to system, and provide no hint as to how a control field can be designed. As such, the possibility for complete 
controllability in a variety of degenerate systems is still unknown and it is widely accepted that coherent control in 
degenerate systems is more difficult than in nondegenerate systems 0, 0, IT^ . Il3| . 

Even more challenging is to establish the complete controllability of a quantum system subject to restrictions on the 
population transfer pathway along which the system may evolve. Shapiro and Brumer [l4j | have shown that there is 
often a loss of controllability if the whole Hilbert space of a quantum system is partitioned and one wishes to generate 
an evolution that docs not pass through some particular Hilbert subspace. Sola et al. [l5j have proposed an optimal 
control theory to maximize the final population in the target state and minimize a time-integral of the population in 
all the unwanted states. Neither of these two previous studies has provided conditions for a class of quantum systems 
subject to strong constraints to be completely controllable. 

Encouraged by rece nt prog ress in extendin g th e stimulated Raman adiabatic passage (STIRAP) 

method and 

its variations 0, 0, ^[ |2fl 0, l23l l24j| for coherent control, we address the controllability of population 
transfer to one of M degenerate states in a system with N > M nondegenerate intermediate states by invoking 
a general adiabatic passage scheme. In particular, we demonstrate that, under certain conditions, it is possible to 
realize complete population transfer from an initial state, via N non-degenerate intermediate states, to one arbitrary 
member of a set of M (M < N) degenerate states, without ever populating any states other than the initial state 
and the target state. The solution is surprisingly simple and general in the sense that realizing such a population 
transfer pathway may not require our knowledge of the exact number of degeneracies or the properties of the M 
degenerate states except for the target state. The solution can also be easily adapted for the creation of an arbitrary 
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superposition of the M degenerate states. The results we have obtained strongly suggest that a class of multi-level 
quantum systems with degenerate states should be completely controllable. The results also imply that complete 
controllability in a class of quantum systems may still be possible even when extremely strong constraints are applied, 
e.g., never populating a Hilbert subspace that is only dimension-two smaller than the whole Hilbert space. 

The central idea of our approach is to manipulate the component fields of the full control field to create an 
(M + N — l)-node null eigenstate of the system dressed by the fields. This multi-node null eigenstate is designed to 
fully correlate with the initial state at early times and then evolve to the target state as the component fields that 
have specified relative phases and amplitudes are turned on and off in a particular order. This scheme is a significant 
extension of our previous adiabatic passage method for the realization of complete control of the population transfer 
branching ratio between two degenerate states [l^. I n addition, both the present study and our previous work 
|19| demonstrate a powerful marriage between weak-field coherent phase control and traditional adiabatic passage 
techniques ^(| that are insensitive to the relative phases of the control fields. 

This paper is organized as follows. In Sec. II we describe our model system. We then present in Sec. Ill our 
adiabatic passage scheme for the complete population transfer from an initial state to one arbitrary member of M 
degenerate states. Simple computational examples that support our theoretical results are described in Sec. IV. We 
discuss our results and conclude this paper in Sec. V. 

II. THE MODEL SYSTEM 

Consider a multi-level quantum system that has a nondegenerate initial state (denoted |0}), and M degenerate 
orthogonal states (denoted I/2}, |/m))- Without loss of generality we denote one arbitrary member of 

the M degenerate states that we wish to transfer population to as |/m)- We assume that the M degenerate states 
cannot be directly coupled with the initial state and that there exist N nondegenerate intermediate states (denoted 
1*2)5 ■ • ■ , Kiv)) that can be coupled with the initial state and some of the M degenerate states. The question 
under consideration is the following: under what conditions can we, in principle, realize 100% population transfer 
from state |0) to state \Jm) without ever populating the intermediate states or any other degenerate states. From a 
mathematical point of view, this question is formidable because (i) the extent of controllability in degenerate systems 
without constraints is not yet resolved and (ii) the extent of controllability in systems subject to strong constraints 
has barely been studied. 




FIG. 1: A schematic diagram of an adiabatic passage scheme for the realization of complete population transfer from an initial 
state |0), via N nondegenerate intermediate states | -ia ) , • • • , |*iv}, to state |/m) that is degenerate with other (M— 1) states 
l/i), I/2), ■ ■ ■ , |/a/-i). Under certain conditions 100% population can be transferred from state |0) to state |/m) without ever 
populating any of the intermediate states or any of the undesired degenerate states. 
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One needs N laser fields to resonantly couple the initial state to the N intermediate states, with the corresponding 
Rabi frequencies denoted by 20 Pl , 2fip 2 , 2ft Pn , and an additional N laser fields to resonantly couple the 
intermediate states with the M degenerate states, with the Rabi frequency associated with state \ik) (k = 1, 2, ■ • • , N) 
and state \fj) (j = 1, 2, • • ■ , M) represented by 2Qs fc . . A schematic diagram of the energy levels and the corresponding 
parameters for the Rabi frequencies is shown in Fig. 1. As in the simplest version of STIRAP ^(|, we assume that 
the laser fields are Gaussian pulses and that they are counter-intuitively ordered. Specifically, the electric field Ep k 
that couples state [0} with state \ik) is given by 

E Pk (t) = E Pk co S (u; Pk t+ cApJ exp [-(t - T) 2 /T 2 ] , (1) 

and the electric field Eg k that couples state with the M degenerate states is given by 

E Sk (t) = E Sk cos(u; Sk t + <f, 8k ) exp [-t 2 /T 2 ] . (2) 

Here t is the time variable, T is the pulse width, E Pk and Eg k are the peak amplitudes of the electric fields, uj Pk and 
Los k are the laser carrier frequencies, 4> Pk and <f>s k are the laser phases. Note that E Pk (t) is delayed by T relative to 
Es k (t). The Rabi frequencies are then given by 

2n Pk =2n Pk exp[-(t-T) 2 /T 2 ], (3) 



2n Skj = 2n Sk] exp [~t 2 /T 2 ] , (4) 
where the peak Rabi frequencies 20, Pk and 2Q,g k . are 

2Cl Pk = ii kE Pk exp(i(j) Pk ), (5) 



2&s kj = fJ>kjE Sk exp(ict>s k ). (6) 

Here fXok = (0|AHfe) an d Mfcj = (*fc|Al/i) are the transition dipole moments. 

In the rotating-wave approximation and in the interaction representation, the Hamiltonian of the system plus the 
component fields of the full control field is given by 
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As suggested by this Hamiltonian, we have assumed that the iV nondegenerate intermediate states will not be coupled 
to one another by the external fields. 



III. ADIABATIC PASSAGE IN DEGENERATE SYSTEMS 



Our previous work on complete control of the population transfer from an initial state to one of two degenerate 
target states 0] suggests that the desired complete population transfer from the initial state |0) to the target state 
|/m) may be achieved by using an eigenstate dressed by the external fields that has nodes on all the states that we do 
not wish to transfer population to. If, by choosing the component fields of the full control field so that a multi-node 
dressed eigenstate is created that correlates first with state |0) and then with the target state |/m) as the fields 
evolve in time, then in the adiabatic limit, i.e., in cases where the dynamics of the system can adiabatically follow 
the evolution of the multi-node dressed eigenstate, this control field will generate complete population transfer from 
state |0) to state |/m)- 
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In particular, we examine here the existence and properties of the null eigenstate (denoted |A)) of H. By definition, 
one has 



H|A> = 0. 



(8) 



In the same representation of H, |A) can be written in terms of its (1 + N + M) components, i.e. 
(z f) ;x 1 ,x 2 ,- ■ ■ ,x N ;y 1 ,y 2 ,- ■ ■ ,y M ) T - Then Eqs. and © give 
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Below we examine the solution to Eqs. ©, (|10|l and llljl in three different cases, specifically, M — N, M < N, and 
M > N. 



A. The M = TV" case 

Let M = N. We denote the M X M matrix S k j (k = 1, 2, • • • , M, j = 1, 2, • • • , M) by S, 

^S 2 i ^s 22 ■ ■ ■ ^S 2M 



Clearly, if 

then Eq. IjlUfl has only one solution 



^Smi ^S M2 ■ ■ ■ ^Smm 

det(S) ^ 0, 
x\ = x% = ■ ■ ■ = xn = 0. 



(12) 



(13) 



(14) 



With this solution for x k (k = 1,2, ••■ ,N), Eq. JJjJ is automatically satisfied. Further, under condition (|13|l . one 
obtains the only solution to Eq. (|llfl . namely 
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where the value of Zq is determined by the normalization requirement and S 1 denotes the inverse of S. Note that 
S _1 can be explicitly written as follows: 
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where A^j is the cof actor of f2s fe . , 



A 



(-l) fc+J 'A fej -, 



(17) 



with being the complementary minor of &s h ■> i-s., the determinant of an (M — 1) x [M — 1) submatrix of S 
obtained by deleting the fcth row and jth column of S. 

It should be pointed out that whether or not det(S) can be zero is a time-independent property inherent to the 
system itself. Using Eqs. |@J and © one obtains 
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Thus, provided that 
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dct(S) will be nonzero at all times. Unless there are some particular system symmetries, it is rare that condition <|19|) 
will be violated. 

Thus, under condition (jl9|) . H has only one null eigenvector with N nodes on the intermediate states and its other 
(M + 1) components zq, y%, y2, ■ • ■ , yu satisfying Eq. I)15|l. Due to the ordering of the laser fields [see Eqs. © 
and Q], this null eigenvector initially correlates with state |0) and then correlates with a superposition of the M 
degenerate states, Ylk=i Vk\fk)- So if the system remains in this null eigenstate, the population will be transferred 
from state |0) to state 2fc=iJ/fe|/fc) without populating any of the N intermediate states as time passes. We must 
still examine whether or not it is possible to realize complete population transfer to one arbitrary memeber of the M 
degenerate states without ever populating the other (M — 1) degenerate states, and if yes, under what conditions. 

Our consideration is based on the relation 
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Hence, if the laser fields are designed such that 
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where is any function of time, then Eqs. I|15|l. H21(l and 1|22J) yield 
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That is, by manipulating the component fields of the total control field we can guarantee that the only null eigenvector 
of H will have additional (M — 1) nodes on states |/i), I/2), • • • , |/m-i}- The total number of nodes of |A) is then 
given by (M + N — 1). It follows that if the dynamics of population transfer adiabatically follows the evolution of 
A), complete population transfer from state |0) to one arbitrary member of the M degenerate states can be achieved 
without ever populating the other (M — 1) degenerate states or the N intermediate states. 
For the specific pulse shape functions given in Eqs. © and condition 1|22() is equivalent to 
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where r\ is a nonzero constant. Using Eqs. JSJ and JSJ), the above condition can be translated to 
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(25) 



Surprisingly, Eq. (|25|) makes no reference to the transition dipole moments that are related to states |/i), I/2}, • • • , 
|/m-i}- Then, even if one has no knowledge about the properties of states I/2), • • • , |/m-i) beforehand, it is 
still possible, by use of adiabatic passage, to completely suppress the population transferred to these states. Note 
also that condition (|25|l requires a definite phase relationship between the component fields of the full control field, 
which is a consequence of the fact that quantum interference effects are directly responsible for the appearance of the 
additional (M — 1) nodes of |A). 

There is another interesting implication of Eq. (|25|l . Let us assume that among the N transition dipole moments 
related to state I of them, say, fJ-2M, • ■ ■ > A*zm, happen to be zero. To still have an (M + N — l)-node 

null eigenvector of H, Eq. [12511 requires Ep 1 = Ep 2 = ■ ■ ■ = Ep l = 0, i.e., removing the I laser fields that couple 
state |0) to the intermediate states |«2), • • • , This clearly results in a simplification of our adiabatic passage 
scheme, but it does not suggest that in these cases the intermediate states |«2), • • • , \ii) are no longer important 
or can be removed from the system, because they are still coupled with the other unwanted degenerate states. These 
intermediate states in this situation can be regarded as a high-dimensional analog of the so-called branch state in a 
five-level extended STIRAP scheme 10], for which case it has been shown that increasing the coupling between an 
unwanted state and the branch state will enhance the yield of the desired state. 



In this case Eq. (fill) requires 
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Under condition (I22II . Eq. 126fl still gives the result of Eq. I|23fl . which can also ensure that Eq. (|27|l is satisfied. It is 
then clear that the null eigenvector, if it exists, will still have (M — 1) nodes on the unwanted degenerate states. 

Next we examine Eqs. © and ijTU|) for M < N. Note first that under condition lf2"2"jl. Eq. © is included in Eq. 
IjlUp. Equation (|10|l indicates that now there are N variables x\, X2, ■ • ■ , a; at and yet the number of the constraints is 
only M. Hence there are more than one solutions to Eq. (|10|l . i.e., H can have multiple null eigenvectors if M < N. 
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The most obvious solution to Eq. i|10|) is given by x\ = x% = ■ ■ ■ = x^ = 0. Under condition i|22|) the associated 
null eigenvector is then given by 

|A) = |Ai> = {z ; 0, 0, • • • , 0; 0, 0, • • • , -z £(i)) T (28) 

This eigenvector is analogous to the only null eigenvector in the M — N case. Let the other null eigenvectors be 
denoted | A2) - Under condition (12211 IA2} can be written as 

|A 2 ) = (z ;xi,x 2> ■ ■ ■ , x N ; 0, 0, • • • ,-z £,(t)) T , (29) 

with X\, X2, • • • ,xn satisfying Eq. I|1U|) and at least one of them being nonzero. Since any linear combination of |Ai) 
and IA2} is still a null eigenvector, we rewrite the other null eigenvectors as 

|A 3 ) = |A 2 )-|Ai) 

= {Q-x u x 2 ,--- ,x N ;Q,Q,--- ,0) T . (30) 

Note that | A3) is orthogonal to Ai. As such, at early times state | Ai) fully correlates with the initial state |0), but 
all the other null eigenvectors represented by | A3) have zero overlap with state |0). 

The important question is whether or not one can avoid populating states | A3) and therefore avoid populating 
the intermediate states during the population transfer. At first glance it seems plausible that since states IA3) are 
degenerate with |Ai), nonadiabatic transitions from |Ai) to | A3) will be unavoidable. Interestingly, this is not the 
case. Specifically, the strength of the nonadiabatic coupling between | A3) and |Ai) is proportional to \(K^\dKi/ dt)\ : 
which turns out to be zero at all times due to the (M + N — l)-node structure of |Ai) and the (M + l)-node structure 
of I A3) [see Eqs. (|28f) and H3U|) ]. Hence, only the null eigenvector |Ai) that has (M + N — 1) nodes is relevant to the 
dynamics of population transfer. 

The analysis here leads to a significant prediction. That is, under condition (|25|l and for M < N, we can always 
achieve adiabatic passage via the null eigenvector |Ai), irrespective of the actual value of M and the values of the 
transition dipole moments between the intermediate states and the unwanted degenerate states. As a result, during 
the entire process of population transfer the undesired degenerate states and all the intermediate states can be never 
populated, even if the exact number of the unwanted degenerate states and their properties are unknown to us. 

Of course, if the exact number of degeneracies M < N is given, it is unnecessary to still use the 2N component 
fields to achieve complete population transfer. As indicated in the preceding subsection, one may remove 2(N — M) 
laser fields and consider only M intermediate states to realize the same control. 

C. The M > N case 

In this case Eqs. @ and IjlUfl in general have only one solution given by Eq. 1141) . So the null eigenvector, if it 
exists, should generically have N nodes on all the intermediate states. However, there are multiple solutions to Eq. 
(|llf) since there are M unknown variables y%, y 2 , ■ ■ • j Um whilst the number of the constraints is only N . Therefore 
H again has multiple null eigenvectors. Denote two orthogonal null eigenvectors of H by |A') and |A"), 

|A') - (4;0,0,--- ,0;y[,y' 2 ,--- ,y' M ) T , 

|A") = (z3,o,...,o ;2/ i',y 2 V--,yM) T (3i) 

The strength (denoted \) of the nonadiabatic coupling between states |A') and |A") is proportional to \(z' )* dz$ / dt + 
^2j = i(y'j)* dy'j I dt\. Due to the normalization requirement, z' Q = y[ = y' 2 = ■ ■ ■ = y' M = is not an acceptable solution 
to Eq. As such, there is no general reason why \ should be zero. In particular, under condition l|22|l . one finds 

that one null eigenvector can be analogous to that in the M < N case, i.e., |A') = (zo; 0, 0, • • • , 0; 0, 0, • • • , — zo£(i)) T . 
Then \ K Iz^dz'^/dt — Zq£* (t)dy'lj /dt\, which is nonzero in general. This being the case, state |A") and therefore 
some of the states |/i), I/2), • • • , |/m-i) will be populated during the population transfer. This is true no matter how 
slowly the laser fields are turned on or off, as states |A") and |A') are degenerate. 

Clearly, then, in the M > N case, our scheme cannot guarantee complete population transfer to state |/m )• However, 
the results above suggest that this problem can be easily fixed if there are more intermediate states available. That 
is, using additional (M — N) intermediate states and adding 2{M — N) laser fields will recover the desired complete 
population transfer. 
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D. Complete population transfer to arbitrary superpositions of M degenerate states 



We note that states I/2), |/m) are just one particular choice of the basis states of an M-dimensional 

degenerate subspace. Indeed, any arbitrary superposition of states I/2}, ■ • • , |/m) can be used as one of the basis 
states of the same degenerate subspace. Consider now a target superposition state 



i/m> = E c ^->' 



(32) 



where Cj (j = 1, 2, • • • , M) are arbitrary coefficients in the superposition state. Let \f' M ) be the last new basis state of 
the same degenerate subspace and let the other new orthogonal basis states be |/{), j/j), • • • , The transition 

dipole moment between state \ik) and state is represented by \J k -. 

We now apply the above general solution to complete population transfer in degenerate systems to state \f' M ) 
instead of \f M ). Then, for M < TV and 
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all population can be transferred from state |0) to the superposition state \f' M ) if the component fields satisfy 
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where 



M 



fi k3 , k = 1,2, • • • ,N. 



(35) 



Since the transition dipole moments fi' k j can be obtained from fikj by considering a unitary transformation between 
two different sets of basis states of the same degenerate subspace, it can be easily proved that condition l|33l) is 
equivalent to condition (|19|l . Note also that there is no need to construct, from superpositions of states |/i), I/2), • ■ • , 
|/m), the explicit forms of states /{), I/2), • • ■ , as their properties are not required to predict that they can 

be never populated during the population transfer. 



IV. NUMERICAL EXAMPLE 



To confirm the validity and feasibility of our general solution to the realization of complete population transfer 
to one arbitrary member of a set of M degenerate states, we present in this section some numerical examples. In 
particular, we consider a system that has N = 7 nondegenerate states as intermediate states and M < N degenerate 
orthogonal states. We define the sum of population in all the intermediate states as P x , the sum of population in all 
the degenerate states except for state \Jm) as P y , and the population in the target state |/V) as P/. For convenience 
all the transition dipole moments are assumed to be real and all the laser phases are set to zero. Hence all the Rabi 
frequencies take real values. To fulfill condition Ij24(l . we assume (lp k = (ls kT , for k = 1, 2, ■ • • , N. 

Figure 2 displays an example of population transfer in the case of M = N. It is seen that during the population 
transfer, the maximum value of P x is less than 0.3%, and the maximum value of P y is less than 0.1%. The final Pf is 
extremely close to 100%. We then arbitrarily alter some of the peak Rabi frequencies that are related to states 
I/2), • • • \fe)- The associated results are shown in Fig. 3. Clearly, population transfer is still almost complete and P x 
and P y again remain negligible at all times. This demonstrates that totally suppressing the population transferred 
to states I/2}, • ■ • , \ f§) does not require us to know their properties beforehand. Similar results are obtained in 
numerous other numerical experiments. 

Figure 4 displays analogous results in a case where two transition dipole moments and [in (and therefore the 
two peak Rabi frequencies &s 17 and r2s 27 ) that are related to the target state I/7) happen to be zero. As predicted 
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FIG. 2: A numerical example of adiabatic passage through a 13-node null eigenstate for the complete population transfer from 
an initial state to one arbitrary member of M = 7 degenerate states. The number of intermediate states is given by N = 7. 
P x represents the population in all the intermediate states, P y represents the population in all the degenerate states except for 
the target state, and Pj represents the population in the target state I/7). The parameters for the peak Rabi frequencies are, 
in units of l/T, given by (fi Pl , Cl P . 2 , &p 3 , fip 4 ,f2p 5 , Q Pe , Cl Pr ) = (60, 90, 60, 120, 90, 99, 135), (f2 Sll , fis 12J &s la , «s 14 , fis 1S! 
Q Sl6 , h Sl7 ) = (90, 15, 0, 150, 36, 18, 60), (Qs 21 , fis 22 , ^s 23 , fis 24I «s 25 , &s 26 , fis 27 ) = (90, 57, 24, 45, 69, 78, 90), (Os 31 , 
£ls 32 , fl S33 , f2 S34; Cl S35 , fl S36 , Oa, T ) = (90, 75, 39, 36, 39, 78, 60), (Cl Sil , fis 42 , n S43 , n S44 , n S45 , n S46 , « S47 ) = (60, 18, 24, 
75, 66, 48, 120), («s 51 , f2 Ss2 , tt Ss3 , fi Ss4 , f2 Ss5 , f2s 56 , Q S57 ) = (39, 27, 93, 15, 66, 78, 90), (fi 5ei) fi Se2 , f2 Se3 , f2 Sfi4 , «s 65 , fis 66 , 
^s 67 ) = (93, 69, 18, 87, 72, 78, 99), and (Q 3n , ps 72 , Clg 73 , n Sri , f2s 75 , Cl Sr6 , «s 77 ) = (36, 54, 48, 57, 96, 78, 135). Note that 
condition 12411 is satisfied since we have chosen Clp k = £ls k7 , for k = 1, 2, • • • ,7. At all times P x < 0.3% and P y < 0.05%. 




time (units of T) 

FIG. 3: Same as in Fig. 2 except fl Sl4 = 39, fis 25 = 39, n S26 = 48, £ls 32 = 45, Q Si3 = 81, f2s 52 = 57, fis 56 = 48, n Sfi2 = 39, 
Qs 72 = 24. At all times P x < 0.2% and P y < 0.03%. 

by Eq. i|25|) . this requires us to remove two laser fields that connect the initial state to two intermediate states. 
Accordingly we set flp 1 and fl p 2 to zero. As seen from Fig. 4, almost complete population transfer with P x and P y 
negligible at all times is also achieved. 

Figure 5 shows a numerical example where all the peak Rabi frequencies that are related to states and I/2) 
are set to zero, while keeping other peak Rabi frequencies the same as those used in Fig. 2. Due to this procedure, 
and I/2) are totally decoupled from the system and effectively there are only five degenerate states. As seen from 
Fig. 5, almost 100% population transfer to state I/7) with P x and P y negligible at all times is also obtained. This 
clearly confirms our previous prediction that our adiabatic passage method works as well even if the exact number of 
degeneracies {e.g., M = 5 or M — 7) is unknown to us. 

To end this section, we note that by increasing the pulse width or the field intensity, or optimizing the time delay 



10 



1 

0.8 



.9 0.6 

H— » 

I 0.4 
0.2 





-2-10 1 2 3 

time (units of T) 

FIG. 4: Same as in Fig. 2 except ft Sl7 = &s 27 = and fi Pl = n P . 2 = 0. At all times P x < 0.03% and P y < 0.01%. 
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FIG. 5: Same as in Fig. 2 except that all the peak Rabi frequencies that are related to states and I/2) are set to zero, i.e., 
Q Shl = &s k2 = 0, for k = 1, 2, • ■ • , 7. At all times P x < 0.3% and P y < 0.05%. 

between the laser fields and the ratio of the peak Rabi frequencies (e.g., Clp k /Ctg kM ), one can always further suppress 
the nonadiabaticity of the dynamics and therefore the population transfer can be made even closer to the adiabatic 
limit. 



V. DISCUSSION AND CONCLUSION 

Controlled population transfer from a selected initial state to a manifold of deg enerate states has previously been 
examined in the context of atomic processes involving magnetic sublevels [21J, l2ll I22L I2IH l24j . Our results are much 
more powerful and general than those previously established. A key difference between this work and previous studies 
is that we have exploited the properties of particular coherent excitations of manifolds with multiple nondegenerate 
intermediate states. In particular, with N nondegenerate intermediate states it becomes possible to adjust, with 2N 
laser fields, all the Rabi frequencies involved in Eq. I|22|l . Hence Eq. <|22l) can always be satisfied by designing the 
laser fields, providing thereby a systematic solution to the design of complete population transfer from an initial state 
to an arbitrary superposition of M < N degenerate states. It is also straightforward to determine if the requirements 
(i.e., condition (|19fl and M < N) are met for our general solution to be applicable. 

Most coherent control approaches based on STIRAP-like dynamics are independent of the relative phases of the 
control laser fields |l7|. As shown in this work and in Refs. and j^, this is not true for all STIRAP-like control 
schemes. On one hand, the phase-dependence of the control process causes even more experimental difficulties; on 
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the other hand the phase-dependence of our adiabatic passage scheme clearly demonstrates that laser phases can be 
also important in STIRAP-likc dynamics and may bring about dramatic opportunities for controlling atomic and 
molecular processes. 

One great advantage of our adiabatic passage scheme is its adaptability in cases where the intermediate states or 
the unwanted degenerate states are not bound states. That is, since the (M + N — l)-node null eigenvector does not 
overlap with any of the intermediate states or any of the unwanted degenerate states, the structure of this multi-node 
eigenvector survives if these states actually decay. As long as the associated decay rate constants are much smaller 
than the peak Rabi frequencies of the component fields |19| , the nonadiabaticity induced by the decaying states will 
be negligible and the dynamics of population transfer would be still given by the evolution of the (M + N — l)-node 
null eigenvector. 

This work extends our previous adiabatic passage method for the realization of complete control of the population 
transfer branching ratio between two degenerate states [l9l | . The extension is from a five-level and four-pulse scheme 
to a (1 + N + M)-level and (2iV)-pulse scheme. While the five- level and four-pulse scheme should be experimentally 
achievable, directly applying the present work to realistic systems would be demanding for M > 3, as it would in 
general require a significant number of laser fields. However, from the point of view of understanding the limits to 
complete controllability in degenerate systems subject to strong constrains and/or with unknown parameters, this 
work is of great theoretical interest. In particular, we have shown that it is possible to realize complete population 
transfer without ever populating a Hilbert subspace that is only dimension-two smaller than the whole Hilbert space, 
even when the exact number of degeneracies and the properties of the unwanted degenerate states are unknown to us. 
We hope that the results of this study will motivate future theretical work on coherent control in degenerate systems. 
For example, it is interesting to consider variations on our adiabatic passage scheme, and to seek alternative, and 
experimentally more feasible, control scenarios that can provide the same type of population transfer pathway that 
is already shown to exist. 

Although our adiabatic passage solution offers a general method for complete population transfer from an initial 
state to an arbitrary member of a set of M degenerate states, there are still two subtle differences between this solution 
and an actual proof of complete controllability of the degenerate system. Firstly, a completely controllable degenerate 
system will have a complete population transfer pathway within a finite time, whereas the complete population 
transfer in our solution, in the most strict sense (i.e., Pf is exactly 100%), can only be realized in the adiabatic limit, 
e.g., for infinitely large pulse width. Secondly, we have used the rotating-wave approximation in constructing our 
solution. Mathematically speaking, such an approximation inevitably modifies the issue of complete controllability. 
It is unclear to what physical extent this approximation changes the controllability of degenerate systems. For the 
above two reasons we regard our adiabatic passage method as a specific physical, but not mathematical, solution to 
the complete control of degenerate systems. 

One potential application of this study is to provide a useful guide for understanding the extremely complex laser 
fields obtained by genetic algorithms in adaptive feedback control of quantum systems 0,0 H^|. Suppose an adapative 
feedback control experiment is carried out in a degenerate system and the control goal is to maximize the population 
in only one of M degenerate states and minimize the population in all the unwanted degenerate states and the 
intermediate states. It would be of great interest to see if the optimized control field suggested by a genetic algorithm 
can be decomposed into 2N > 2M laser fields with their frequencies, amplitudes, and relative phases close to what 
is suggested by our adiabatic passage scheme. If this is the case, then such a feedback control experiment becomes 
a realization of our general solution to complete population transfer in degenerate systems and more applications of 
this work can be expected. 

To summarize, we have shown that under certain conditions complete population transfer from an initial state, via 
N nondegenerate intermediate states, to an arbitrary superposition of M (M < N) degenerate states is achievable 
by adiabatically following an (M + N — l)-node null eigenstate that is created by designing the relative phases and 
amplitudes of the component fields of a control field. The results may find applications in quantum information 
processing in atoms and molecules and shed considerable light on the issue of complete controllability in degenerate 
systems under strong constraints and/or with unknown parameters. 
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